Conformal invariance, multifractality, and finite-size scaling at Anderson localization 

transitions in two dimensions 
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We generalize universal relations between the multifractal exponent ao for the scaling of the typical 
wave function magnitude at a (Anderson) localization-delocalization transition in two dimensions 
and the corresponding critical finite size scaling (FSS) amplitude A c of the typical localization length 
in quasi-one-dimensional (Q1D) geometry: (i) When open boundary conditions are imposed in the 
transverse direction of Q1D samples (strip geometry), we show that the corresponding critical FSS 
amplitude A° is universally related to the boundary multifractal exponent ag f° r the typical wave 
function amplitude along a straight boundary (surface), (ii) We further propose a generalization of 
these universal relations to those symmetry classes whose density of states vanishes at the transition, 
(iii) We verify our generalized relations [Eqs. (^) and (Q)] numerically for the following four types of 
two-dimensional Anderson transitions: (a) the metal-to- (ordinary insulator) transition in the spin- 
orbit (symplectic) symmetry class, (b) the metal-to- (Za topological insulator) transition which is 
also in the spin-orbit (symplectic) class, (c) the integer quantum Hall plateau transition, and (d) 
the spin quantum Hall plateau transition. 

PACS numbers: 73.20.Fz, 05.45.Df, 72.15.Rn 
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I. INTRODUCTION 

Localization-delocalization (LD) or Anderson localiza- 
tion transitions of non-interacting electrojas-,are contin- 
uous phase transitions driven by disorder au When dis- 
order is weak, the single-electron wave functions are ex- 
tended over the whole sample. Sufficiently strong disor- 
der localizes electrons within a finite region in space. The 
linear size of this region is the localization length £ char- 
acterizing the typical size of the wave functions ijj(r)xl As 
the disorder strength is reduced, the localization length 
increases and eventually diverges at an LD transition 
point. The localization length is the analogue of the cor- 
relation length at non-random continuous phase transi- 
tions. At the LD transition point, wave function- ampli- 
tudes obey scale-invariant, multifractal statistics;Btl that 
is, the disorder-averaged g-th moment of the square of the 
absolute value of wave function has a power-law depen- 
dence on the linear dimension L of the system, with an 
exponent that is a non-linear function of q.Qrl 

Let us recall that continuous phase transitions in non- 
random systems are known to be quite generally de- 
scribed by conformally-invariant field theories. Confor- 
mal symmetry is especially powerful in two dimensions 
(2D), where its presence leads to an infinite number of 
symmetry constraints. This, in many cases, allows fpr. 
a rather complete description of critical properties. EjEII 
Effective (field) theories describing the random LD tran- 
sitions are also expected to possess conformal symmetry. 
In fact, we have recently shown by numerical simulations 
of a standard LD transition occurring in two dimensions, 
namely of the metal-insulator transition in the 2D spin- 
orbit (symplectic) symmetry class,t3that multifractal ex- 



ponents of critical wave functions evaluated on a straight 
boundary and those at a corner are related thtojugh a 
simple relation dictated by conformal symmetry.! 3 ! 

Conformal symmetry is known to impose strong con- 
straints on finite-size scaling (FSS) for phase transi- 
tions in non-random systems with quasi-one-dimcnsional 
(Q1D) geometry. For these systems Cardy has shownlifl 
that the correlation length £ which characterizes the de- 
cay of the tyn-ppint correlation function of any (confor- 
mal primaryEJiij) operator along a cylinder or a strip of 
width M, is related to the bulk or surface (x s ) scal- 
ing dimension of the operator in two dimensions through 
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2irxb, cylinder (periodic BC), 
ttx s , strip (open BC). 



(1) 



Here BC stands for boundary conditions imposed in the 
transverse direction. 

The generalization of Eq. (Q) to scale-invariant disor- 
dered 2D systems was first provided in the study of ran- 
dom 2D diluted ferromagnets in Ref. |l(| [for the 2D bulk 
exponents and Q1D cylinder geometry (periodic BCs)]. 
In a random system the scaling of an observable (such 
as, for example, a "spin") is in general characterized by 
the set of scaling dimensions x q of all its i%th moment 
disorder averages. Equation ([j]) generalizes^ to all these 
moments. In particular, the correlation length £ 9 char- 
acterizing the exponential decay of the q-th moment of 
a correlation function of the observable in Q1D cylinder 
geometry is related to the 2D scaling exponent by 

— = 2irx qi cylinder (periodic BC). (2) 

At the same time, by using an expansion about q = 



2 



of the g-th moments in the 2D system and in the Q1D 
cylinder geometry, it was demonstrated in Ref. [fl6| that 
such a relationship holds also for the corresponding u typ- 
icaV quantities referring to a fixed disorder realization. 
In particular, if ao denotes the typical^ 2D bulk scaling 
dimension of the observable, and if l/£ denotes the Lya- 
punov exponent characterizing the inverse of the typical 
Q1D correlation length in cylinder geometry, then again 
the relationship 
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cylinder (periodic BC) 



(3) 



holds. 

Later, Refs. f§,|l7j proposed a corresponding formula 
in the context of LD transitions in two dimensions, 



M , h 
T = Mao 



2), 



(4) 



[the shift by two between the r.h.s. of Eq. (Q) and of Eq. 
(|J) arises from different conventions]. Here, £ p is the typ- 
ical Q1D localization length in cylinder geometry (the 
subscript p of £ p denotes periodic BCs imposed in the 
transverse direction). The exponent a\ in Eq. (^J) char- 
acterizes the scaling of a typical critical wave function 
amplitude in the bulk of a 2D system of linear dimension 

R, 



In 



-a b n \nR, 



(5) 



where the overbar stands for the disorder average. Equa- 
tion (|]) has been confirmed numerical ly for , the integer 
quantum Hall (IQH) plateau transitionotZlta and for the 
2D metal-insulator transition in the spin-orbit (symplec- 
tic) symmetry class. E3E§EJ 

We note that the relation ([|) , in the form presented, is 
only valid for systems in which the average bulk density of 
states (DOS) is constant and non-vanishing at the tran- 
sition. This is the case for LD transitions in the three 
Wigncr-Dyson classes. These include the IQH plateau 
transition and the LD transition in the spin-orbit (sym- 
plectic) class. However, as is now well known, there are 
symmetry classes in which the DOS vanishes at the tran- 
sition. This is the case, for example, for the so-called spin 
quantum Hall transition of the Bogoliubafedc. Gennes 
(BdG) quasiparticlcs in symmetry class Cc2rE3 (in the 
nomenclature of Ref. |2^) . 

In this paper we derive a generalization of the relation- 
ship ([|) between the exponent «g and the typical Q1D 
correlation length £ p for LD transitions in 2D with a van- 
ishing critical DOS. The result is 



— =2ir(a b -2 + x p ), 



(6) 



where the exponent x p characterizes the critical behavior 
of the (bulk) DOS (x p = in the Wigncr-Dyson classes). 

Furthermore, we derive a FSS formula for the typical 
Q1D localization length, when open BCs are imposed 



in the transverse direction. The specific open BC we 
consider in this paper is a reflecting BC which means 
that the system simply ends at the boundary, so that 
there is no current flowing across the boundary. The 
second line of Eq. (Q) suggests that the localization length 
should be related to a surface exponent characterizing 
multifractality of critical wave functions near boundaries 
of disordered systems.Ej Indeed, our result is the formula 
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n(a s -2 + Xp), 



(7) 



where now a§ is the surface (i.e., boundary) exponent 
characterizing scaling of a typical wave function ampli- 
tude near a straight (reflecting) boundary, is defined 
in the same way as ajj in Eq. (^|), except that now the 
point r is close to a straight boundary of the 2D sys- 
tem of linear dimension R. The typical QID localization 
length £ D is computed in the geometry of a strip of width 
M with open (reflecting) BCs imposed in the transverse 
direction (the subscript o stands for "open"). 

The organization of this paper is as follows. In Sec. II 
we derive Eqs. (^|) and (J?]). In Sec. Ill we verify both these 
equations numerically by computing the critical FSS am- 
plitude (A p or A ) of the typical QID localization length, 
defined as 



p M ' 



M 



(8) 



for both types of BCs (the factor 2 in this definition 
is standard convention). We verify Eq. (R) for (a) the 
metal-to- (ordinary) insulator transition in the spin-orbit 
(symplectic) class [class All of Rcf. , (b) the LD tran- 
sition between a metal and a Z2 topological insulator in 
the "quantum spin Hall" (QSH) effected which also be- 
longs to the spin-orbit (symplectic) class [class All of Rcf. 
|25| , (c) the IQH plateau transition in the unitary symme- 
try class [class A of Ref. . [The bulk relation, Eq. (|Fj) , 
was already verifkyl .farj-siystcms (a)-(c), where x p = 0, 
in previous work.BEJOcj] We finally verify numerically 
Eqs. (|) and (g) for the spin quantum Hall transition in 
symmetry class C of Ref. |25| Table | summarizes the nu- 
merical results presented in detail in Section III. Section 
IV presents our conclusions. 



II. LOCALIZATION LENGTH AND 
MULTIFRACTALITY 

In this section wc provide a derivation of Eqs. (||) and 
(0) . Let us begin with a brief discussion of the underlying 
assumptions. We are interested in scaling properties of 
the disorder average of some physical observable [e.g., the 
local DOS (LDOS)] at an LD transition point. One can 
recast this disorder average into a statistical average of a 
properly defined operator O in a certain field theory (e^™ 
a replica or supersymmetric nonlinear sigma model) Jala 
The scaling properties of O at the critical point are then 
controlled by the fixed point of the rcnormalization group 
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TABLE I: A list of ao, x p and A c for various universality classes. The values of ao marked by * are from the references listed 
in the last column. Those without * are obtained in this paper. The fifth column shows A c calculated from Oq using Eq. (^) 
and from using Eq. (R), combined with Eq. These values of A c should be compared with those obtained from fitting 
(as explained in Section III) and shown in the sixth column. 

system BCs ao x p A c from ao A c from fit Ref. 

symplectic (M-I) open 2.429 ± 0.006 1.48 ± 0.02 1.50 ± 0.01 this paper 

symplectic (M-QSH) open 2.091 ± 0.002* 7.00 ±0.15 7.20 ± 0.01 § 

IQH open 2.385 ± 0.003 1.654 ± 0.013 1.624 ± 0.002 this paper 

SQH in class C periodic 2.137* 1/4 0.8225 0.8189 ± 0.0004 [Es| 

SQH in class C open 2.326* 1/4 1.105 1.101 ±0.002 \24\ 



(RG) flow of the corresponding field theory. We are now 
ready to state the-lwo important assumptions we make 
in our dcrivation:EZI 

• The fixed-point theory is a conformal field theory. 

• At the fixed point of the RG transformation, the 
operator O is a primarylI3 field operator in the con- 
formal field theory. 



A. Finite-size scaling in cylinder geometry, and 
bulk exponents 

Let us consider a disordered electronic system at its 
critical point, confined to a disk of radius R in the 2D x-y 
plane, or equivalently, the complex plane with the coor- 
dinate z — x ± iy. We assume that all along the bound- 
ary of the disk there is a metallic electrode attached, 
thus allowing for the electron in the system to escape. c3 
This (absorbing) boundary condition introduces a finite 
broadening r\ of the single-particle levels in the system. 
We assume that the broadening is of the order of the 
mean level spacing in the system. This provides a regu- 
larization for Green's functions and the LDOS as follows: 



G ± {z,z';E) = Y r ^ z) ^ Z '\ 



Pe(z) = — [G+(z,z;E)-G.(z,z;E)} 



-£hM*)l 

IT * » 



(E - E n ) 



(9) 



(10) 



(z) of the closed system are 
\ilj(z)\ 2 d 2 z = 1. The in- 



Here the wave functions ip r , 
normalized in the disk: f, . 

tegral of the LDOS pe{z) over the disk gives the global 
DOS pe multiplied by the disk area nR 2 . 

Statistical properties of metallic or critical wave func- 
tions at energy E are closely related to those of the 
LDOS.Ej In particular, if we are interested in the scaling 
of the moments of such wave functions and the moments 
of the LDOS, we can write symbolically 



Pe{z) 
nR 2 p E 



(11) 



Disorder averages of powers of the LDOS pe(z) (as well 
as those of products of Green's functions) are represented 
by expectation values of operators in the corresponding 
field theory.Elta We denote this by 



(12) 



where the angular brackets denote the expectation value 
in the field theory. Here O q (z) is the operator which 
corresponds to the q-th moment of pe{z). (We point out 
that here and in what follows the power q can take any 
real values.c3) In view of Eq. ([n]), the same operator 
represents moments of the wave function iPe(z)' 



[R 2 PE ) q \^ E (z)\ 2 ^ (O q (z)). 



(13) 



Notice that the global DOS is self-averaging and can be 
pulled out of the disorder average along with powers of 
the radius R. The product R 2 pe oc <5 _1 , where 5 is the 
mean level spacing in the disk. 

Now we concentrate on the wave functions and the 
DOS at the critical energy, E = E c , and drop the sub- 
script E. The global DOS p may vanish at criticality 
in the infinite system. In a finite system the disorder- 
averaged p always has a power-law behavior 



R- 



(14) 



where the exponent x p vanishes in the standard Wigner- 
Dyson classes but is known to be non-zero in other sym- 
metry classes. For example, at thC|-(2.D.) spin quantum 
Hall transition in symmetry class C,clrE3 the exact value 



mmc 



is known: x p = 1/4 
We now ■— .make 
assumptionsE3 of canformal 



of 



the previously stated 
invariance and the fact 
that O q is a primarytLEI conformal scaling operator with 
the bulk scaling dimension x h q at the LD transition. If 
we choose a point \z\ <S R close to the origin of the disk, 
then the one-point function (the field theory expectation 
value) scales as 



(O q (z)) ~ R~ 



(15) 



Combining this with Eqs. (13) and (|14|), we obtain the 
scaling of the moments of the critical wave functions: 



\ip(z)\ 2 i - R 



-2q-x +qx p 



(16) 
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for \z\ <C R. Notice that the exponent of R on the right 
hand side should vanish at q = 0, and should be —2 at 
q = 1 due to the normalization of the wave function. 
These conditions determine 



x 



0, 



x l 



(17) 



Some important details of the definition and properties 
of multifractaLjexponents are in order here. A slightly 
more detailedEII ( "coarse-grained" ) description of multi- 
fractal wave functions (in 2D) involves breaking the sys- 
tem into little square boxes Bi of size r x r labeled by 
i. The number of these boxes N scales as N ~ (R/r) 2 . 
One then calculates the probability pi for an electron to 
be in the i-th box as 



Pi 



m z )\ 2 d 2 z, 



(18) 



B, 



and forms the so-called average generalized inverse par- 
ticipation ratios 



N 



(19) 



(We have assumed that the system is homogeneous af- 
ter disorder average.) Equation ([!(]) implies the scaling 
relation 



2(q-l)+x b q -qx p , (20) 



where the set of exponents T q is usually referred to as the 
multifractal spectrum. 

Note that the probabilities pi whose moments enter 
the definition of P q are bounded by ^ p, ^ 1. This 
bound implies that P q must be a non-increasing func- 
tion of q, since pf 1 ^ pf for q\ < q 2 . Moreover, since 
p q = cxp((7lnp) is convex as a function of q, the same is 
true for P q . Then the multifractal spectrum r q in Eq. (^p| ) 
must be a non-decreasing, concave function of q. Gener- 
ally speaking, there may be a value of q = qj where r q 
has a horizontal tangent. Then it follows that r q = const 
for q ^ qj. Such change in the behavior of r q from an 
increasing function to a constant is often referred to as 
"freezing" or "termination" (see Ref. || for more details) . 
In all known cases the value qj where such termination 



occurs satisfies g/ > 0. Then we can safely use Eq. (|16|), 
and similar equations in the following sections, in the 
vicinity of q = without worrying about a possible ter- 
mination transition. 

Expanding both sides of Eq. (|l^) in q about q = 
yields the typical scaling exponent, Eq. (H), where 



dx b q 
dq 



q=0 



Next, let us consider the conformal mapping 

'2tt 



w 



M 

— In z, 
2tt 



exp 



(21) 



(22) 



which maps the disk to the semi-infinite cylinder of cir- 
cumference M in the complex w-plane, 



w = u + iv, u ^ L 



M 
2^ 



lni?, 0^v<M, (23) 



with an absorbing boundary condition at u = L. The. 
assumption that O q is a primary conformal operators 
allows us to relate its expectation value on the cylinder 
to that in the disk: 



(o q H) 



dz 



dw 
(-) 

Km) 



q (O q (z)) 
exp 



2- 
M 



(24) 



This immediately gives the moments of the LDOS in the 
cylinder: 



IpM}" 



exp 



—x b JL-u) 
M qy ' 



(25) 



From the exponential decayiiU of the moment [p(w)] 9 
away from the end of the semi-infinite cylinder in Eq. 
(p5|), for sufficiently small positive values of q 7 we iden- 
tify the "g-dependent localization length" £ p (q) in the 
cylinder geometry as 



M 
2nx b q 



(26) 



(Here 'p' denotes again the 'periodic' BCs of the cylin- 
der.) The typical Q1D localization length £ p in cylinder 
geometry is read off from the typical exponential decay 
of the LDOS away from the end of the semi-infinite cylin- 
der: 



In p(w) = — 



\L-u\ 



Expanding again Eq. (E5J) in q about q = yields 



M 



2tt- 



dx b 
dq 



q=0 



= 27r(ao - 2 



(27) 



(28) 



where we have used Eq. (|2l]). This is our previously 
mentioned result, Eq. (g), which generalizes Eq. (Q) to 
all symmetry classes, including those with critical DOS. 

In Section HID we numerically verify Eq. (^) for the 
spin quantum Hal effect (symmetry class C) by comput- 
ing numerically the FSS amplitude A p = 2£ p /M of the 
typical Q1D localization length £ p in cylinder geometry; 
according to our above-obtained result (p8[) this quantity 
is predicted to equal 



A„ = 



1 



P Tt{a b Q -2 + x P Y 



(29) 



with x p = 1/4. 



5 



B. Finite-size scaling in strip geometry, and 
surface (boundary) multifractal exponents 

We now apply the same arguments to discuss finite- 
size scaling in the presence of open (reflecting) BCs in 
the transverse direction (strip geometry). 

For this purpose we first consider the operator O q 
placed close to the origin in the interior of the half disk 
\z\ ^ R, Imz 0. The boundary of the system on the 
real axis is assumed reflecting, and the rest is attached to 
a metallic lead, as in the previous section. In this situa- 
tion, the expectation value of O q (z) for \z\ <C R is given 
by@ 



(O q (z)) ~ R- 



(30) 



where the boundary scaling dimension x s q (the super- 
script s stands for "surface") is typically different from 
the bulk dimension x q . In analogy with Eq. (|l6| ) we now 
have, upon making again use of Eq. (|l3|), 



\lp(z)\ 2 1 ~ R- 2 «-**+V*r, 



(31) 



where the same exponent x p (a bulk exponent) enters 
through the global DOS. Note that Eq. @ still implies 
= 0, but now, in the boundary case, there is no re- 
striction on x\ (in contrast to the bulk case: see Eq. ( |l6| ) 
and the subsequent text). Also, in complete analogy to 
the bulk case, the exponent of R in Eq. (|3l]) must be a 
monotonic function of q. Upon expanding both sides of 
Eq. ( |3l"| ) in q about q = 0, one obtains the scaling ex- 
ponent c>!g of the typical wave function amplitude at the 
boundary, 



ln|V(*)| 5 



where now 



2 x n 



dx q 
dq 



q=0 



(32) 



(33) 



Next, in order to relate this to the strip geometry, we 
use the conformal transformation 



M 

w = — In z, 

7T 



exp 



KM J 



(34) 



which maps the half disk to a semi-infinite strip of width 
M in the w-plane: 

M , 
w = u + iv, u^L = — In R, < d < M. (35) 

7T 

The expectation value_pn the strip now follows again 
since O qi as a primaryEj conformal operator, transforms 
simply under conformal transformations, 



(o g H) 



dz 



dw 



'(O q (z)) 



(-Y 9 

VaiJ 



exp 



M< {L - U) 



(36) 



From this we obtain the exponential decay of the mo- 
ments of the LDOS away from one end of the strip, 



[p{w)] q 



exp 



u) 



(37) 



As in the bulk case, the exponential decay0 of the right 
hand side in Eq. (|37|), for sufficiently small positive values 
of q, gives the "g-dependent Q1D localization length" 
along the strip 



t ( \ M 
irx q 



(38) 



As before, the typical Q1D localization length £ Q in strip 
geometry is obtained by expanding both sides of Eq 
in q about q = 0, 



M 
To 



dx s q 
dq 



q=0 



ir(a s Q 



(39) 



where we have used Eq. (B3). This is our previously- 
announced result from Eq. (p). 

In subsequent sections we verify Eq. (0) for various LD 
transitions by computing numerically the FSS amplitude 
A Q = 2£ Q /M of the QID typical correlation length £ D on 
the strip ('o'='open', reflecting BCs) which, according to 
our result, is predicted to be equal to 



A, 



?r(ag - 2 + x p ) 



(40) 



III. NUMERICAL RESULTS 

In this section we present the results of our numerical 
simulations supporting Eqs. (|^) and (|4Cj|). For conve- 
nience, we have gathered all the relevant fittingparame- 
ters and other numerical data in a single table |J. 

In this section we have to distinguish off-critical and 
critical values of the QID localization lengths, £ and £ c , 
and the corresponding FSS amplitudes, A and A c (for 
both periodic and open BCs). All £ and A that have 
appeared in the previous sections denoted values at the 
critical point. 



A. Spin-orbit (symplectic) symmetry class 

To compute the localization length at the LD transi- 
tion in the syjsiplectic class, we employed the so-called 
SU(2) model,L3 a tight-binding model on the square 
lattice, with random on-site disorder and fully random 
SU(2) hopping. 



1. Localization length (strip geometry) 

We obtained the typical localization length from the 
smallest Lyapunov exponent of transfer matrices for very 



G 



TABLE II: A list of parameters obtained or used in the FSS analysis for the scaling functions defined in Eqs. (^l|) and (^?]). 
Here A c , v, and y are obtained through fitting. Nd and N p denote the numbers of data points and fitting parameters used in 
the fitting procedure, respectively. The fitting functions are truncated at the orders P and Q. \ 2 an d g denote the the values 
of chi squared and the goodness of fit probability, respectively. 



system 



BCs scaling function 



A r 



N d N p P Q x 2 



<1 



symplectic (M-I) reflecting Eq 

symplectic (M-QSH) reflecting Eq 

IQH reflecting Eq 

SQH in class C periodic Eq 

SQH in class C reflecting Eq 



(«) 
(41) 

(«) 
(57) 
(il) 



1.50 ±0.01 2.79 ±0.03 -1.03 ±0.03 85 9 2 2 86.2 0.2 

7.20 ±0.01 - -0.81 ± 0.08 8 3 7.2 0.2 

1.624 ±0.002 2.55 ±0.01 -1.29 ± 0.04 134 6 3 2 144.0 0.2 

0.8189 ± 0.0004 1.335 ±0.016 -0.94 ± 0.01 73 8 2 2 56.2 0.7 

1.101 ±0.002 1.335 ±0.005 -1.05 ±0.02 93 9 3 2 86.1 0.4 




long Q1D lattices. We imposed hard-wall, i.e., reflect- 
ing BCs in the transverse direction and hence our Q1D 
samples had strip geometry. Our systems had a maxi- 
mum size M = 128 in the transverse direction. Figure 
|l](a) shows the FSS amplitude A D = 2£ Q /M of the typical 
Q1D localization length as a function of the on-site disor- 
der strength W for various system sizes M and at fixed 
energy E — (band center). The curves for the var- 
ious system sizes intersect at different points reflecting 
largcJinitc-size effects, in contrast to the case of periodic 
BCsH 

To determine the critical value of the FSS amplitude 
A o c , we performed a FSS analysis incorporating correc- 
tions to scaliag arising from the leading irrelevant scal- 
ing variablc.E-3 Specifically, we took a scaling function for 
the FSS amplitude of the form A = F( X M 1 / v XM y ), 
where \ is the relevant scaling variable, and £ is the 
leading irrelevant scaling variable whose scaling expo- 
nent y < 0. The exponent v characterizes the diver- 
gence of the 2D localization length £ upon approaching 
the LD transition point, £ ~ x" ■ We expanded the scal- 
ing function around the critical point W = W c , setting 

x = (w- w c )/w c , 



p=l 



f Q 

K = K.c ±E«p (x Ml/ "Y + MV E b « (* Ml/ " 

9=0 



(41) 

We fitted the numerical data to Eq. (|4lj) with P = Q = 2 
by taking W c , a p , b q , is, and y as fitting parameters. We 
obtained 



W r . 



(42) 



6.192 ±0.007, A 0jC = 1.50 ±0.01, 
v = 2.79 ±0.03, y = -1.03 ±0.03. 

The details of the fitting are summarized in Table 0. 

These results are in good agreement with those ob- 
tained by Asada et al. for the SU(2) modelta with peri- 
odic BCs: W c = 6.199 ± 0.003 and v = 2.75 ± 0.04 at 
E = 0. The good quality of the fit can be seen from the 
scaling collapse, shown in Fig. 0(b), of the data for the 
corrected FSS amplitude A Q defined by 



A n = A n 



M y b. 



q (xM 1 



(43) 
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FIG. 1: (Color online) (a) Dependence of A on W at E = 
for various values of M. A at E = 0. (b) Scaling plot of 
A at E — - see Eq. (jiil). The obtained parameters are 
A , c = 1.50 ± 0.01, W c = 6.192 ± 0.007, v = 2.790 ± 0.025, 
y = -1.026 ± 0.03, ai = -1.69 ± 0.03, a 2 = 0.70 ± 0.02, 
bo = 1-24 ± 0.03, bi = -2.36 ± 0.08, and b 2 = 4.7 ± 0.3. 



2. Surface multifractal exponent ao 



In our previous publicationEj we reported the value 
ckq = 2.417 ± 0.002 for the surface exponent, which was 
obtained from numerical simulations on L x L lattices of 
system sizes up to L = 120. We performed averaging over 
more than 6 x 10 4 disorder realizations. The lattices had 
periodic BC imposed in one of the two directions, but 
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open BC in the other direction, so our system had the 
geometry of a finite cylinder. Here wepupdate the value 
for a.Q reported in our previous workJI3 We use larger 
system sizes up to L = 180, and average over up to 10 5 
disorder realizations. 

The surface exponent was obtained from the system 
size dependence of the wave function amplitude in the 
vicinity of the boundary, according to 



« In |^(a:)| 2 » ~ -oglnL + c,. 



(44) 



Here x = <D(L ), I > 1, and c is a constant of order 
L° . The double angular brackets represent both ensemble 
average and spatial average along the boundary of the 
cylinder in each disorder realization. First we tried a 
linear fitting to Ecl_(44|) of our numerical data for the left 
hand side of Eq. (44j7usmg system sizes 24 ^ L ^ 180, 
with two fitting parameters and c. This resulted in 
the value 



= 2.4195 ±0.0013. 



(45) 



Substitution of this value into Eq. (^0() gave 

A , c = 1.518 ± 0.005. (46) 

This analysis, however, ignored corrections from irrele- 
vant scaling variables and was not quite correct, since we 
now know from the previous subsection that such correc- 
tions are appreciable for the FSS amplitude A for open 
BC. We therefore re-analyzed the data, assuming scaling, 
with corrections from the leading irrelevant variable.Ej 
We define 



A(x) := - 



((ln|^)| 2 )) 



1 



InL 



where we take y = —1, as suggested by Eq. (42). The 
fitting of the same data to Eq. (|47l) yielded 



which leads to 



A , c 



2.429 ±0.006, 



1.48 ±0.02 



(48) 



(49) 



with the help of Eq. (fig). We see that the A D:C obtained 



from the transfer matrix method (42) is consistent with 



these results. The value of reported in Eq. ( |48| ) has 
larger error bars, which needs to be improved in future 
numerical work. 



B. Metal to Z2 topological insulator transition in 
quantum spin Hall systems 

The Z2 topological insulator is a time-reversal invari- 
ant topological insulator in two dimensions, which pos- 
sesses a topologically protecteeLKramers pair of extended 
edge states at its boundaries.E3 The Z 2 topological insu- 
lating states can be realized in materials with strong spin- 
orbit interactions, as evidenced by recent experiments 



on HgTc/(Hg,Cd)Te quantum wclls.Ea In the presence of 
disorder, this system undergoes a two-dimensional metal- 
insulator transition from a Z2 topological insulator to a 
metal, as one changes the Fermi energy. On symme- 
try grounds, this LD transition is expected-to belong to 
the spin-orbit (symplectic) symmetry class.EH Indeed, the 
critical exponent v for the diverging localization length 
(a bulk property) at the metal to I2 topological insulator 
transition is found_.to agree with the value obtained for 
the SU(2) model,E3 which describes the metal to (ordi- 
nary) insulator transition in this symmetry class. Similar 
agreement is found for the multifractal exponents for crit- 
ical wave functions in the bulk.ES However, the multifrac- 
tal exponents characterizing wave function amplitudes at 
the sample boundary turn out to be different at the two 
metal-insulator transitions. 

Here we show that, at the metal to Z2 topological in- 
sulator transition, the FSS amplitude A o c (Eq. (||)) for 
the typical Q1D correlation length in strip geometry, is 
related by conformal invariance to the boundary multi- 
fractal exponent at the same transition. 



1. Localization length (strip geometry) 

To compute the localization length at the metal to Z2 
topological insulator transiti on, w e employed the quan- 
tum spin Hall network model.EZIfia An important parame- 
ter in this network model is the one controlling the prob- 
ability of tunneling at the nodes of the network, which we 
denote by X. The numerical results shown below were 
obtained at the critical point X c = 0.971 with &41y ran- 
dom SU(2) spin rotation symmetry on each link.E3 Figure 
U shows the dependence of the FSS amplitude A Q (M) := 
2£ (M)/M of the typical Q1D localization length £ Q (M) 
on a strip of width M (M = 8, 10, 12, 16, 24, 32, 48, 64). 
Here M is the number of nodes of the network model in 
the transverse direction across the Q1D strip. This cor- 
responds to transfer matrices of size AM x AM. In order 
to find the critical value A 0jC of the FSS amplitude A„ 
in the large M limit, we assumed that A Q at X = X c 
has a power-law finite-size correction due to a leading 
irrelevant variable with dimension y < 0: 

A (A = X c ) = A OjC + & M y - (50) 

Fitting the data to this form (see Fig. ||), we obtained 

A 0:C = 7.20 ±0.01 (51) 

with y = -0.81 ± 0.08 and b = -1.0 ± 0.1. The details 
of the fitting are summarized in Table |l|. 

2. Surface multifractal exponent ctQ 

The surface multifractal exponent at the metal to Z2 
topological insulator transition was obtained in Rcf . [p8| . 
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0.05 0.1 



II M 

FIG. 2: M dependence of A at the metal to topological 
quantum spin Hall insulator transition. The solid curve is a 
fit to Eq. @ with A D , C = 7.20 ± 0.01, y = -0.81 ± 0.08, and 
b = -1.0 ±0.1. 

By using larger system sizes this value was recently im- 
proved in Ref. Q to 

a s Q = 2.091 ±0.002. (52) 

Substituting the improved value into Eq. ( [io|) yields the 
FSS amplitude 

K . c = 7.00 ±0.15. (53) 

This value is consistent with Eq. (|5l|). The larger error 
bar in Eq. (B3h results from the fact that the denominator 
in Eq. @ (with x p = 0) contains a s - 2 = 0.091 ± 
0.002. Neither of the numerical analyses in Refs. |38|, p9[ , 
used to obtain Eq. (|52]), included effects of the leading 
irrelevant variable, in contrast to Eq. (|5l]). These effects 
may influence the value of cvq and possibly result in better 
agreement with Eq. (|5l"l). 

C. Plateau transition in the integer quantum Hall 
effect 

To compute the localization length £ Q and the surface 
multifractal exponent ctp at the plateau transition in the 
IQH effect, |-a5e-. employed the Chalker-Coddington net- 
work modelc30 in strip geometry with M nodes in the 
transverse direction across the strip. This corresponds to 
transfer matrices of size 2M x 2M. The plateau transi- 
tion is reached by tuning a parameter 9 which controls 
the tunneling probability at the nodes of the network 
model. For this model the critical value 9 C is known ex- 
actly. 

1. Localization length (strip geometry) 

The typical localization length £ D in Q1D strip geom- 
etry was computed numerically from the smallest Lya- 
punov exponent of the transfer matrices. The largest 
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FIG. 3: (Color online) (a) Dependence of A on the node 
parameter 9 in the Chalker-Coddington model of the strip 
geometry. The vertical dashed line indicates the critical point 
6 = 6 C . Inset: M dependence of A at 6 = 6 C ; the solid curve is 
a fit to Eq. (JHoj) . (b) Scaling plot from FSS analysis including 
corrections from the leading irrelevant scaling variable. The 
parameters used for the plot are v — 2.55 ±0.01, ai = 2.518 ± 
0.016, a 2 = 2.179 ±0.027, a 3 = 1.393 ±0.051, b = 1.26 ±0.7, 
bi = 2.016 ± 0.086, and b 2 = -0.73 ± 0.26. 

system size (the number of network model nodes in the 
transverse direction) that we studied, was M — 64. Fig- 
ure ||(a) shows the FSS amplitude A Q = 2£ /M of the 
typical localization length as a function of the network 
model tunneling parameter 9 for various transverse sys- 
tem sizes A/.EJ For 9-> 9 C the network model is in the 
quantum Hall phasc.l!3 As seen from Fig. ||, the crossing 
point of the curves moves towards 9 = 9 C as M increases, 
indicating the presence of finite-size corrections. To find 
the critical value of the FSS amplitude A D of the typical 
Q1D correlation length in the large M limit, we fitted 
the data to Eq. ( [so)) [see the inset of Fig. 0(a)], to obtain 

A , c = 1.624 ±0.002, (54) 

y = -1.29 ± 0.04, and b = 1.26 ± 0.7. The details 
of the fitting are summarized in Table [n} Figure ||(b) 
shows the data collapse from the FSS analysis using Eqs. 
( J4l| ) and (El) with \ = (# - d c )/0 c and the values of 
A 0jC , y, and 6o obtained above. This FSS analysis also 
yielded v = 2.55 ± 0.01 for the critical exponent of the 
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FIG. 4: The Chalker-Coddingt on netw ork model on a cylin- 
der. In the notation of section III C 1 L = 3 on this figure. 
There are 4L 2 = 36 links, and the unitary evolution operator 
U is a 36 x 36 matrix. 



diverging (2D bulk) localization length, which is close 
to the value obtained recent large-scale numerical 



study, v = 2.593 ± 0.006 



2. Surface multifractal exponent a§ 

The surface multifractal exponent ccq at the plateau, 
transition was recently obtained by the pxesent authorsE3 
and by Evers, Mildenberger, and Mir lino It was found in 
these works that the multifractal analysis for the Chalker- 
Coddington model suffers from large finite-size correc- 
tions. To reduce these corrections, we have used, in the 
multifractal scaling analysis in Ref. [Q, numerical data 
obtained only for large system sizes. Here we used an al- 
ternative approach by taking into account corrections to 
scaling arising from a leading irrelevant scaling variable 
using Eq. @. 

The geometry of the Chalker-Coddington network 
model that we used is shown in Fig. I There are two 
types of nodes forming two sub-lattices (denoted A and 
B in the figure), such that the A sublattice has the size 
L x L (L = 3 in the figure). The links of the network 
form zigzag shaped rows and columns; there are 2L such 
rows and 2L such columns so that the total number of 
links is 4L 2 . Integer x and y coordinates are assigned 
to the centers of links. We imposed periodic BC in the 
vertical y direction, and reflecting BC in the horizontal x 
direction. The links in the first and the last columns at 
x = 1 and x = 2L are called the edge links. The discrete 
time evolution of wave functions defined on links of the 
network model is governed by a unitary evolution oper- 
ator U for one discrete time step, which is determined 
by the [Scattering S matrices at the nodes of the network 
modelo In our case this operator is a 4L 2 x 4L 2 uni- 
tary matrix. For each disorder realization, we obtained 
one critical wave function that is the eigenvector of U at 
9 = 9 C and whose eigenvalue is closest to unity among 
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l/ln(2L/0, l/ln(2M) 

FIG. 5: (Color online) (a) Spatial dependence of logarithm of 
probability density ((In \ip(x)\ 2 )). The system size is changed 
as L = 15, 20, 25, 30, 35, 40, 50, 60, 80, 120, 180 from the top to 
the bottom, (b) Dependence of A(x) := — ((In \ip(x)\ 2 )}/ InL 
on the effective system size L/l at x = 1 (squares, I — 1), 
2 (open circles, I — 1) and at edge plaquettes with coarse- 
graining (filled circles, 1 — 2), where I is the (linear) size of 
boxes used to define the coarse-grained wave function proba- 
bilities. The solid curves are the fits to Eq. (M). Also shown 
by crosses is the dependence of 2[1 + 1/ttA\1T — 9 C )] on the 
width M and its fit (dashed curve) to Eq. (p0|). 



all the eigenvectors. The largest system size we studied 
was L = 180, and the disorder average was taken over 
3 x 10 5 realizations for L ^ 60, over 5 x 10 5 realizations 
for L = 80, and over 2 x 10 5 realizations for L = 120, 180. 

Figure ||(a) shows the x dependence of ((In |?/>(x)| 2 )), 
where the double angular brackets stand for both the 
average over disorder realizations and the spatial aver- 
age along the periodic y direction. We clearly observe 
in Fig. ||(a) Friedel-like oscillations near the edges of the 
cylinder, which become less pronounced as L is increased. 
(Such oscillations are absent in the SU(2) model dis- 
cussed in the previous section.) Figure |^(b) shows how 
A(x) = —((\n\ip(x)\ 2 }}/lnL approaches a constant value 
with increasing L at the left boundary (x = 1,2). The 
solid curves show the fitting of A(x) to Eq. (|47]) at x = 1 
(squares) and x ~ 2 (open circles). To minimize the 
corrections coming from the Friedel-like oscillations, we 
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defined the coarse-grained wave function amplitude on 
each plaquette and calculated the corresponding A for 
the plaquettes along the edge (shown as red filled cir- 
cles). Fitting this coarse-g rained data to Eq. ( fl7|) with 
y = —1.29 obtained in Sec. Ill C 1 yielded 



a% = 2.385 ±0.003, 



(55) 



where the error bars reflect only statistical errors. This 
result is consistent with that of Ref. Qj (ag = 2.386 ± 
0.004). Figure |(b) shows that fitting of A(x = 1) and 
A(x = 2) gives similar values of ag. Substituting Eq. ( |55| ) 
into Eq. (|40|) yields 



A, 



1.654 ±0.013, 



(56) 



which should be compared with K o c = 1.624 ±0.002 [Eq. 
(|S4|)] obtained from the transfer matrix calculation. As 
we see in Fig. |5](b), finite-size corrections to A and A Q 
are still quite large at L = 180. This makes the extrapo- 
lation of these quantities to L — > 00 difficult; we cannot 
exclude the possibility of having systematic errors in ad- 
dition to the statistical errors included in Eqs. ([34]) and 
( p6| ) . Given the presence of this uncertainty, we conclude 
that our numerical results are consistent with Eqs. (Q) 
and <M). 



D. Spin quantum Hall plateau transition of BdG 
quasiparticles in symmetry class C 

In this section we discuss the verification of Eqs. @ 
and (^) for symmetry class C, which is known to possess 
a vanishing critical DOS (x p > 0). In our simulations 
we used an appropriate generalization of the Chalker- 
Coddington network modelJi3 which we refer to as the 
class C network model. This model has a control pa- 
rameter e (in the notation of Ref. [^|), and is critical 
at e = 0. Exact values for critical exponents, v = 4/3 
and x p = 1/4, w^iA obtained through mapping toxlassi- 
cal perpDlationnd'E3 The exact values of the bulkEa and 
surfaceEj multifractal wave function exponents x^ s are 
also known at q = 2,3. However, exact results for the 
FSS amplitudes of the typical Q1D correlation lengths, 
A PjC and A o c , and the typical wave function scaling ex- 



ponents a ' s are not available. 



1. Localization length (cylinder and strip geometries) 

We numerically obtained the FSS amplitudes of the 
typical Q1D localization length of the class C network 
model for both cylinder and strip geometries. A previ- 
ous numerical studytD of FSS of the typical localization 
length in cylinder geometry did not report the value of 
A PjC . Here we present results for the FSS amplitudes A PiC 
and A 0jC corresponding to cylinder and strip geometries, 
respectively. 
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FIG. 6: (Color online) (a) Dependence of A p on e for several 
values of M in the class C network model of cylinder geometry. 
A critical point is known to be located at e c = 0. (b) Scaling 
plot of A p , obtained after subtracting corrections to scaling 
from a leading irrelevant scaling variable. The parameters 
used for the plot are A p , c = 0.8189±0.0004, v = 1.335±0.016, 
y = -0.94 ± 0.01, a 2 = -1.66 ± 0.10, a 4 = 3.64 ± 0.33, 
60 = 0.185 ± 0.003, and b 2 = 0.58 ± 0.31. (c) M dependence 
of A p at e — 0. The solid curve is a fit to Eq. (|5c|). 



Cylinder Geometry: Figure ^(a) shows the dependence 
of the FSS amplitude A p of the typical Q1D correlation 
length on the parameter e for various values of the trans- 
verse width M, obtained in cylinder geometry. The FSS 
amplitude A p is symmetric about the critical point e c = 
when periodic BCs are imposed. Hence in the FSS anal- 
ysis we have to use an expansion in even powers of e, 



A — A 



p 
P =i 



2p 



Q 

9=0 



2q 



The result of fitting of the data in Fig. ||(a) to Eq. ( |57| ) 
is shown in Fig. |6|(b). The dependence of A p (e c ) on the 
width M at the critical point e c = is plotted in Fig. 
|6|(c). We obtained 



A PiC = 0.8189 ±0.0004 



(58) 



and v = 1.335 ± 0.016. The details of the fitting are 
summarized in Table 0. The latter result is consistent 
with the exact value v = 4/3, indicating good accuracy 
of our numerical results. 
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FIG. 7: (Color online) (a) Dependence of A on e for several 
values of M in the class C network model of strip geometry. 
The critical point is located at e c = 0. (b) Scaling plot of A 
from the FSS analysis with subtraction of corrections from 
a leading irrelevant scaling variable. The parameters used 
for the plot are A , c = 1.101 ± 0.002, v = 1.335 ± 0.005, 
y = -1.05 ± 0.02, ai = 2.225 ± 0.024, a 2 = 3.221 ± 0.083, 
a 3 = 2.91 ± 0.20 6 = 0.960 ± 0.015, bi = 1.846 ± 0.095, and 
b 2 = 2.26 ± 0.54. (c) The M dependence of A D at e = 0. The 
solid curve is a fit to Eq. 



Strip Geometry: Figure ^(a) shows the FSS amplitude 
of the typical Q1D correlation length in strip geometry. 
With reflecting BCs imposed in the transverse direction, 
the model possesses edge states for e > (thC|-spin quan- 
tum Hall phase, possessing topological order) Jl3 Since A„ 
is not a symmetric function of e, we use the FSS function 
in Eq. (E5) . Figures 0(b) and (c) show the result of the 
FSS analysis and the width M dependence of the FSS 
amplitude A Q (e c ) for the typical correlation length in the 
strip, respectively. From this analysis we obtained 



A, 



1.101 ± 0.002, 



(59) 



and v = 1.335 ± 0.005. The details of the fitting are 
summarized in Table [TI]. 

2. Multifractal exponent ao 

The bulk and surface multifractal exponents cvg and 
an for the class C network model have been obtained 



a b ~ 2.137, 



a s ~ 2.326. 



(60) 



Substitution of these values into Eqs. ( p9| ) and 
spectivcly, with x p = 1/4 yields 

A PiC = 0.8225, A . c = 1.105. 



(61) 



These values are consistent with the values presented in 
Eqs. (|58|) and (pi3) obtained by our FSS analysis. 



IV. CONCLUSIONS 

In this paper we have generalized the formula relating 
the multifractal exponent ao of the typical wave function 
amplitude in a 2D sample to the FSS amplitude A c of the 
typical localization length in a Q1D sample. Our gener- 
alization is twofold, resulting in Eqs. (^) and (Q). Our 
Eq. ([|) extends the relation to unconventional symmetry 
classes where the global density of states vanishes at crit- 
icality. Our Eq. (0) extends the relation to the case when 
the Q1D sample has strip geometry, instead of cylinder 
geometry which was always considered in earlier stud- 
ies. In this case the multifractal exponent describes 
the scaling of typical wave function amplitude near the 
sample boundary. 

We have verified generalized Eqs. (|^) and (|7]) numeri- 
cally for systems in four different universality classes: (a) 
the metal-to-insulator transition in the spin-orbit (sym- 
plectic) symmetry class, (b) the metal-to-(Z2 topological 
insulator) transition also in the spin-orbit (symplectic) 
class, (c) the integer quantum Hall plateau transition, 
and (d) the spin quantum Hall plateau transition. Our 
numerical results are summarized in Tables [| and ||. 
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It is knownu that the exponents x h q are positive only within 
some range < q < q, . Outside this range x b q < 0, which 
leads to LDOS moments [p(u>)] 9 that are exponentially 
growing away from the cylinder boundary. This can be ex- 
plained by the following argument. In cylinder geometry, 
the wave functions are localized, and the level broadening 
rj created by the metallic lead attached at the edge u — L 
is strongly coordinate dependent. Indeed, the overlap of a 
wave function localized at a point w away from the cylin- 
der edge with the extended wave functions in the lead is 
exponentially small, which results in an escape rate, that 
is, a level broadening that is exponentially decaying away 
from the edge: n{u) <x e~( L ~">/£. Here £ is the typical 
localization length. Thus, away from the edge, the profile 
of the LDOS (as a function of energy) consists of well re- 
solved narrow Lorentzian peaks (see Eq. (lid)) of width of 
order n{u) and height proportional to l/n(u). Only within 
the range < q < g* corresponding to x\ > are the 
moments of the LDOS determined by the minima of the 
LDOS profile (where the LDOS is proportional to ri(u)). 
For other values of q the dominant contribution to the mo- 
ments of LDOS comes from the maxima in the LDOS pro- 
file. In the described situation the behavior of the moments 
of the wave functions is very different from that of the mo- 
ments of the LDOS. However, we can directly relate the 
latter to the transmission through the cylinder from the 
metallic lead. Essentially, the transmission will be expo- 
nentially small away from the edge of the cylinder if the 
energy of the incoming wave (in the lead) corresponds to a 
minimum in the LDOS profile, but will be of order one (it 
cannot be bigger) if the energy is in resonance with one of 
the broadened energy levels (that is, a peak in the LDOS 
profile). 

The dependence on z, which encodes the distance from the 
boundary, is suppressed here. 
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In this phase A D increases with M, due to the presence of 
edge states, which increase the conductance (similar to the 
case of a metallic phase). 
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